Excitation spectrum of a toroidal spin-1 Bose-Einstein condensate 
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We calculate analytically the excitation spectrum of a toroidal spin-1 Bose-Einstein condensate 
that is subjected to a homogeneous magnetic field and contains vortices with arbitrary winding 
numbers in the m_F = ±1 components of the hyperfine spin. The spectrum can be tuned by 
varying external parameters, such as the strength of the magnetic field. A rotonlike spectrum can 
be obtained, or an initially stable condensate can be made unstable by adjusting the magnitude 
of the magnetic field or trapping frequencies. The structure of the instabilities can be analyzed by 
measuring the particle density of the spin components. 
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Bose-Einstein condensates (BECs) confined in toroidal 
traps have been subject to many experimental studies 
recently This research cover topics such as the 

observation of persistent current [E, phase slips across 
a stationary barrier stochastic [3| and deterministic 
Q phase slips between vortex states, the use of toroidal 
condensates in interferometry [H[ , and the stability of su- 
perfluid flow in a spinor condensate [(§. These experi- 
ments have given rise to theoretical studies discussing, 
e.g., the excitation spectrum and critical velocity of a 
superfluid BEC Q and the simulation of the experi- 
ment using the Gross-Pitaevskii equation 0, Q and 
the truncated Wigner approximation Most of the 
experimental and theoretical studies concentrate on the 
properties of persistent currents. The phase of a toroidal 
BEC changes by 2nk as the toroid is encircled, the in- 
teger k being the winding number of the vortex. In a 
singly-connected geometry a vortex with > 1 is typi- 
cally unstable against splitting into vortices with smaller 
fc. In a multiply connected geometry this process is sup- 
pressed for energetic reasons. In Ref. Q it was shown 
experimentally that a vortex with winding number three 
can persist in a toroidal single-component BEC for up 
to a minute. In other words, toroidal geometry makes 
it possible to avoid the fast vortex splitting taking place 
in a singly-connected BEC and study the properties of 
vortices with large winding number. 

In this letter, we calculate the linear excitation spec- 
trum of a toroidal quasi one-dimensional (ID) spin-1 
BEC. Motivated by the experimental results of Refs. 
[!, @ , we assume that the splitting of vortices occurs on 
a very long timescale in a spinor condensate where only 
one spin component is populated. The dominant insta- 
bilities can then be assumed to arise from the spin-spin 
interaction. For related theoretical studies on toroidal 
two-component condensates, see, for example, [Io|, 11 1- 
In our analysis the population of the mp = component 
is taken to be zero initially, making it possible to calculate 
the excitation spectrum analytically. This type of state 
is straightforward to prepare experimentally and makes 



it possible to observe the proliferation of instabilities by 
measuring the density of the mp = component. 

The order parameter of a spin-1 Bose-Einstein conden- 
sate is ijj = (ipi, ipo, ip-i) T , where T denotes the trans- 
pose, and it fulfills the identity ip*ip = n^p>, where n^D 
is the total particle density. The condensate is con- 
fined in a toroidal trap given in cylindrical coordinates 
as U(r, z, (p) = m [uj 2 (R — r) 2 + uj 2 z 2 ] /2, where R is the 
radius of the torus and <jJ r ,to z are the trapping frequen- 
cies in the radial and axial directions, respectively. We 
assume that the condensate is quasi ID, so that the or- 
der parameter factors as ip(r, z, ip; t) = vb r . z (r, z)ip v (ip; t), 
where ip r - z is complex-valued and time independent. If 
the system is exposed to a homogeneous magnetic field 
oriented along the z-axis the energy functional becomes 



dip (-e-^ - - pF z + qF 2 ^j ^ 



(1) 



where e = h 2 /2mR 2 , fj. is an effective one-dimensional 
chemical potential, and F = (F x ,F y ,F z ) is the (dimcn- 
sionlcss) spin operator of a spin-1 particle. The magnetic 
field introduces the linear and quadratic Zeeman terms, 
given by p and q, respectively. The sign of q can be 
controlled experimentally by using a linearly polarized 
microwave field 12j. The strength of the atom-atom in- 
teraction is characterized by go = 4irh 2 (ao + 2a,2) /3m and 
cj2 = A:Txh 2 {a2 — ao)/3m, where clf is the s-wave scatter- 
ing length for two atoms colliding with total angular mo- 
mentum F. The scattering lengths of 87 Rb used here are 
ciq = 101. Sas and 02 = 100. 4as fl3j | . measured in units 
of the Bohr radius as- For 23 Na the corresponding values 
are ag = 50. Oa^ and a-i = 55. las [Til ]. The magnetiza- 
tion in the z-direction, f z = J dip4 ! l,( l p)F z tp lp (ip)/N, is 
a conserved quantity; the corresponding Lagrange mul- 
tiplier can be included into p. Here N is the condensate 
particle number. We assume that in the initial state the 
spin is parallel to the z axis. In [151 ] it was argued that 
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in a homogeneous system the most unstable states are 
almost always of this form. The initial state reads 



1 + /, 



,0,e i9 e ik - iv 



1-fz 



(2) 



where 9 is the relative phase and the integer k±\ is the 
winding number of the nip = ±1 component. The en- 
ergy and stability are independent of 9 and therefore 
we set 9 = 0. If fc x = 1 and = 0, Vll describes 
a half-quantum vortex (Alice string), see, e.g., Refs. 

The populations of i/ju are time independent 
and the Hamiltonian reads H\\ = (gon — fj,) I + (g 2 n f z — 
p e «)F z + q e «F 2 , where p cS = p - e(kf - k 2 _ 1 )/2 and 
?cff = q + + k 2 _ 1 )/2. The time evolution operator 
of ipu is U\\(t) = e ~ ltH u/ h . We calculate the excitation 
spectrum in a basis where ip\\ is stationary (IBI . [l9j . In 
this basis, the energy of an arbitrary state ip is given 

by E?™[iP] ee E m [U\itp] + ih(il>\ V) ^> and thc 
time evolution of ip can be obtained from ihdip/dt = 
5E^[iP}/8^. We write ip((p;t) = ip\\(tp) + Sip(tp;t) and 
expand the time evolution equation to first order in Sip. 
The perturbation is written as 

oo 

Si> s (!p;t) ee e^Y, u rAt)e lsv -v*, s (t)e- lsv , (3) 

s=0 

where j = 0, ±1 and fco = 0. Straightforward calculation 
gives the equation 



ih 



a 



w 



dt V w i;-i 



Bz®Bi 



w 

Wlj-l 



(4) 



where © indicates the direct sum of matrices and 
the subscript i gives the dimension of Bi. Only Bi 
depends on time. The vectors wo and Wi ; _i are 
defined as w = (u 0;s , Uo ;s -fci-fc_i ) T and Wi ; _i = 
(ui-g, U-i-g, Vi- S , V-i- tS ) T . The matrix elements of B2 read 

(-82)11 = es 2 +g 2 n, 

(-82)22 = e(s - ki - /s_i) 2 + g 2 n, 

and (^2)21 = — (-62)12- The operator B4 can be written 

as 



B 4 = 



h + x 

X 



D 



-X 
-es 2 l 2 - X 



D 



(8) 



where I2 is the 2x2 identity matrix, D 
2es diag(fci, fc_i), and X is defined as 



X llt2 2=2g+(l±f,), X 12 = 2g-^l-fl (9) 

and X 2 i = X 12 . Here g± ee (g ± g 2 )n/A and + (-) 
refers to Xu (^22). The 6x6 matrix B 2 © -B4 is in 



a block diagonal form, the blocks being the B 2 and B4 
matrices. It is then possible to calculate the eigenvalues 
and eigenvectors of B 2 and -B 4 independently. B 2 affects 
only the mp = component, that is, the magnetization 
remains unchanged but the direction of the spin vector is 
modified. The corresponding modes are therefore called 
spin modes. For similar reasons, the eigenvalues of B4 
are called magnetization modes. 

The eigenvalues of B4 can be calculated straightfor- 
wardly but they are too long to be shown here. The 
eigenvalues are independent of q and can be written as 
Hu)j(s) = 2esfc+ + Hu)j(s), where uij depends on fc_ = 
(fci — fc_i)/2 but is independent of fc+ ee (fci + fc_i)/2. 
Consequently, modes with differing fc+ but equal fc_ have 
identical stability. We denote the real and imaginary 
parts of ujj by and wj, respectively Magnetization 
modes are unstable if at least one ujj is positive. If f z = 
the eigenvalues simplify and are given by 



^1,2,3,4(5) \ ft=0 = 2es/c H 



(10) 



±2* Us 2 



-{s 2 +4fc 2 ) +g + ± ^ek 2 _(es 2 +4g+) + g 2 _ 



The signs are defined such that ++, — h, +— , and 

correspond to 0Ji,ui 2 ,uj3, and W4, respectively. Unstable 
modes appear when the term inside the square brackets 
becomes negative. For rubidium and sodium g + > 0, 
guaranteeing that lo\ and u 2 are real. Only W3 can have 
a positive imaginary part, see Fig. [T] The amplitude 
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FIG. 1. The amplitudes of the unstable spin and magneti- 
zation modes for rubidium and sodium. Here e — Q.75\g2\n, 
q = 0, and the unit of lo^ 5 is \g2\n/h. The lines have been 
drawn by treating s as a continuous parameter; dots indi- 
cate the actual allowed non- vanishing values of w^. In (c) 
and (d) the curves are reflection symmetric with respect to 
s = k+ = (fci + fc_i)/2. 

of the unstable modes grows as \k-\ increases. Now the 
allowed values of s are integers; the modes correspond- 
ing to s — are always stable, but unstable modes are 
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present for s = 1,2,..., \_J 4Ar?_ — 2g2n/e\ . Here 
indicates the floor function. It is therefore not possible 
to choose the parameters in such a way that a single 
mode s > 1 is unstable. A lower bound for the value of 
e yielding at least one unstable mode in a sodium BEC 
is given by the equation e(4fc^ — 1) > 2g2n. Because 
e > 0, the magnetization modes of a sodium conden- 
sate (g2 > 0) with fc_ = and |fc_| = 1/2 are always 
stable. This is visualized in Fig. 0Jb), where w 3 (s) cor- 
responding to (fci,fc_i) = (0,0) and = (2,1) 
is seen to vanish for every s. In a rubidium condensate 
(.92 < 0) with fc_ = unstable modes exist if e < 2\g2\n; 
if \k-\ > 0, instabilities are present regardless of the 
value of e. For both rubidium and sodium the wavenum- 
ber of the fastest-growing instability is approximately 

/2A 



given by the integer closest to ^/2/3^J 4k^_ — 2g2n/e. As- 
sume next that there is one dominant unstable mode. 
This is in particular the case if only the s = 1 mode 
is unstable. Using the eigenvector corresponding to 
u; 3 (s) |/, = o w e nn d that the position-dependent popula- 
tion difference AF z (ip;t) = pi(ip;t) — p_i(<p;t), where 

Pi = \bP\\)j + <%| 2 ! is g iven b Y 



AF z (<p; t) = ae"3* cos(6>) + /3e 2 ^* sin(20) 



(11) 



where = arg[ui ;s (0)]/2 + s<p + uj' 3 t and a and j3 are 
constants that depend on s, fc_, and the interaction coef- 
ficients. The presence of unstable magnetization modes 
can be experimentally confirmed by measuring AF z ((p). 
If fc_ =0 we obtain f3 = 0, leading to AF Z having s 
peaks as a function of the coordinate if. For nonzero /c_ 
both a and /3 are in general nonvanishing and AF Z is a 
sum of two competing terms with periodicity s and 2s. 
If o; 3 (s) > we get w 3 (s) = 2esk + , so that the nodes 
of AF Z rotate around the torus as time evolves if k + is 
nonzero. We see that although the value of k + does not 
affect the stability, it can alter the behavior of the density 
perturbations resulting from the instabilities. 

We now turn to the spin modes. The time dependence 
of B2 can be eliminated by a simple change of basis. In 
the new basis, the eigenvalues giving the spin modes read 



f"^5,6(s) = 2ek + (s — k + ) 



(12) 



± yj{e[(s - fc+) 2 - k 2 _ ] + g 2 n - q}~ - (1 - f 2 )(g 2 n) 2 , 

where + (— ) corresponds to uj§ (loq). If k + = 0, the 
effect of vortices can be taken into account by scaling 
q — > q c g = q + efc^, i.e., the spin modes of a system 
with (fcx, fc-i) = (A;, —k) and q = q are equal to the spin 
modes of a vortex- free condensate with q = q + ek 2 . Spin 
modes are unstable if and only if the term inside the 
square root is negative. Now only W5 can have a positive 
imaginary part. The fastest-growing unstable mode is 
obtained at e[(s — fc+) 2 — fc 2 ] + g2n — q = and has the 
amplitude ?ko\{s) = |g2|ri\/l — f z - Unlike in the case 



of the magnetization modes, the maximal amplitude is 
bounded from above and is independent of the winding 
numbers. By adjusting the strength of the magnetic field, 
the fastest-growing unstable mode can be chosen to be 
located at a specific value of s, showing that it is easy to 
adjust the stability properties experimentally. At f z = 
the width of the region on the s-axis giving positive W5 



is |y k 2 _ + q/e — J k 2 _ + q/e — 2g2n/e\. This region can 
thus be made narrower by increasing e, or q. Since 
the magnetization modes are insensitive to the magnetic 
field, the properties of the spin and magnetization modes 
can be tuned independently. The winding number depen- 
dence of unstable spin modes is illustrated in Figs. QJc)- 
(d) . Interestingly, by tuning e and q a rotonlike spectrum 
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FIG. 2. The real (1^5) and imaginary (015) component of ojs 
(spin mode) for (a) rubidium and (b) sodium. Here k\ = 
— fe_i, k 1 is an arbitrary integer, and e = 0.3|c?2 |n.. For the 
blue solid and blue dotted lines q c fi = q + ekL = 5.5|(?2|h and 
for the orange dashed line q c R — — 2\g2\n. The unit of uf^ 1 is 
\g2\n/h. The blue solid and orange dashed lines show uj\ and 
the blue dotted line gives 1^5 . These lines have been drawn by 
treating s as a continuous parameter; dots (circles) indicate 
the actual allowed non- vanishing values of uf 5 (u£,). 

can be realized, see the solid and dotted blue lines in Fig. 
[5J Now the phonon part of the spectrum is missing, but 
the roton-maxon feature is present. For f z = k + = 
the roton spectrum exists if q > max{0, 2g2n}. Because 
only integer values of s are allowed, it may happen that 
the imaginary parts W3 and W5 are nonzero only in some 
interval of the s axis that does not contain integers, see 

Figs. DIc)-nid) an( l l^S- EI a )- l n thi s case rotonic exci- 
tations are stable. Alternatively, there can be unstable 
modes close to the roton minimum, see Fig. [2Kb) and 
Ref. [2(J. As evidenced by the orange dashed lines in 
Fig. [21 the roton spectrum can be made to vanish simply 
by decreasing q. It is known that a rotonlike spectrum 
can exist in various types of BECs, such as in a dipolar 
condensate (see, e.g., [2ll423| ). in a Rydberg-excited con- 
densate or in a spin-1 sodium condensate prepared 
in a specific state (2^. In the present case the rotonlike 
spectrum exists both in a sodium and rubidium BEC and 
the state [Eq. ([2])] giving rise to it is easy to prepare ex- 
perimentally. Note that the roton-maxon feature exists 
also in a vortex- free condensate and for any \f z \ < 1. 
These results suggest that the roton-maxon character of 
the spectrum is rather a rule than an exception in spinor 
BECs. By changing the value of q an unstable spin mode 
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can be introduced or made to vanish. For example, us- 
ing the parameter values corresponding to the blue solid 
line in Fig. [^a), we find that by decreasing (increas- 
ing) the value of q e g from 5.5\g2\n to 4|<? 2 |" (6|g2|")> the 
s = 4 (s = 5) modes can be made unstable. This opens 
the way for quench experiments of the type described in 



25l |26| . In the present system the number of excited un- 
stable modes and their wavenumbers can be controlled 
with a high precision. Instead of altering q, instabilities 
can also be induced by making e smaller by changing the 
trapping frequencies. 

The properties of unstable spin modes can be studied 
experimentally by measuring the density and phase of 
the mp — spin component. Assuming that there is one 
dominant unstable spin mode located at wavenumber s, 
we find the density of the mp = component to be 



\u\ 2 + \v\ 2 



Po { V ;t) = \8Mv;t)\ 2 ~^ t {\i 

- 2\uv\ cos[6 u + 9 V + 2u\t + 2{s - k+)ip]\, (13) 

where u = u - s (0), v = Uo;a-fei-fe_i(0), 9 U is defined 
through u = \u\e lSu , and 9 V is defined similarly. The 
number of nodes in po is 2|s — k+\, that is, if 2k+ is even 
(odd), only an even (odd) number of nodes can exist. 
The density peaks resulting from an instability rotate 
around the torus if w§ = 2efc+(2s — fc + ) is nonzero. The 
ratio of the components of the eigenvector Wo = (u, v) T 
corresponding to W5 reads v/u = (a — 6)/yl — f^gin, 
where a = e[(s — k + ) 2 — k_^\ + g^n — q and b = 
yV -(1- f 2 )(g 2 n) 2 . If there is a dominant unstable 
mode with wavenumber s we get \v/u\ = 1 and Eq. (|13p 
gives po(ip;t) ss 4|u| 2 e 2 "5 t sm 2 [u!lt+ (s — k + )ip], where we 
have set 9 U = 9 V = for simplicity. The density van- 
ishes at (p = (cjjjt + mt)/(k + — s), where n is an arbitrary 
integer. From the fact that 5ipo ~ sin[o;§t + (s — h+)<p] 
we see that the sign of Stpo changes at every point where 
the density vanishes. This holds true regardless of the 
value of f z . The ±7r phase jump at the density minima 
resembles the behavior of the phase of a dark soliton [27| ■ 
We look next at the values of e corresponding to two 
recent experiments. Using the Thomas-Fermi (TF) ap- 
proximation and the parameter values of the sodium 
experiment Q we get e ss 0.0452"- This means that 
magnetization modes are stable if |fc_| < 3, but spin 
modes are unstable even when \k—\ = 0. In the follow- 
ing k- is such that only spin modes are unstable. If 
< q < 0.0452", fz = 0, and = (0,0), the insta- 

bility appears as a position-independent, homogeneous, 
increase in pq. On the other hand, if (k±, fc_i) = (1, 0) we 
get po(ip;t) ~ e 2a; =* sin 2 [(e£ + <p)/2], so that the density po 
is peaked around tp = —et+ir. This illustrates that even a 
small e can lead to a strongy winding number dependent 
density pq. The parameters of the rubidium experiment 
Q give e w 0.18|<72|"- Now some magnetization modes 
are unstable regardless of the values of winding numbers. 



If q e B ~ 2I52I" the s = 4 spin mode, situated close to 
the minimum of the roton spectrum, is the only unstable 
spin mode. At q c g — 0, the s — 1,2,3 modes are unsta- 
ble, and at q e g = — 2\g 2 \n all spin modes are stable. Using 
the TF approximation we get e ~ {Lo r Lu z N R i )~ 1 / 2 , show- 
ing that the stability of the magnetization modes can be 
controlled by adjusting the trapping frequencies, number 
of particles, and the radius of the toroid. 

The first experimental realization of a toroidal spin-1 
BEC was reported very recently @. The stability of a 
rubidium BEC with vortices with equal winding numbers 
in the tuf = 1 and mj? = components was found to 
depend strongly on the population difference of the two 
components, the most unstable situation corresponding 
to equal population. Although not directly comparable, 
our analysis agrees qualitatively with this result: The 
amplitude of unstable spin and magnetization modes in- 
creases as f z — > 0. 

In conclusion, we have calculated analytically the lin- 
ear excitation spectrum of a toroidal spin-1 BEC. In 
this calculation the strength of the magnetic field and 
the winding numbers of the nip = ±1 components were 
treated as free parameters. We found that the shape 
of the spectrum can be tuned by adjusting the strength 
of the magnetic field. In particular, a spectrum with 
a roton-maxon structure can be realized both in rubid- 
ium and sodium condensates. Furthermore, by varying 
the magnetic field, an initially stable condensate can be 
made unstable. We have also shown that some unstable 
modes lead to a dark soliton -like wavefunction of the 
mp = spin component. Finally, by considering two re- 
cent experiments on toroidal single-component BECs, we 
argued that various types of instabilities can be realized 
in a spin-1 condensate having the same parameter values 
as these experiments. 
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